Relativistic mean field study of the properties of Z=117 nucleus and the decay chains of 293 ,294 117 

isotopes 
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We have calculated the binding energy, root-mean-square radius and quadrupole deformation parameter for 
the recently synthesized superheavy element Z=117, using the axially deformed relativistic mean field (RMF) 
model. The calculation is extended to various isotopes of Z=l 17 element, strarting from A=286 till A=310. We 
predict almost spherical structures in the ground state for almost all the isotopes. A shape transition appears 
at about A=292 from prolate to a oblate shape structures of Z=117 nucleus in our mean field approach. The 



most stable isotope (largest binding energy per nucleon) is found to be the 
a-decay Q a and the half-lives T t 
magic numbers at N=172 and/ or 184 



1 17 nucleus. Also, the Q-value of 
are calculated for the a-decay chains of 293 117 and 294 117, supporting the 



PACS numbers: 21.10.Dr, 21.10.Ft, 21.10.Gv, 21.10.Tg 
I. INTRODUCTION 

Nuclei can survive beyond the macroscopic limit, far into 
the transuranium region, where the necessary balance between 
the nuclear and Coulomb force is achieved only through shell 
stabilisation. Superheavy element (SHE) are hypothesised to 
exist in this region. The next double shell closer, beyond 
208 Pb, predicted at Z=l 14, N=184, may have suprisingly long 
half-life, even of the order of a million year |Q]-|5[]. 

Experimentally, till to-date, elements upto Z= 1 1 8 have been 
synthesised by heavy ion reactions (a [ZD, with half-lives 
ranging from a few minutes to about a milli-second. The 
more microscopic theoretical calculations have predicted the 
next region of stability, beyond Z=82, N=126, as Z=120, 
N=172 or 184 (§{11 and Z=124 or 126, N=184 JH 0. 
However, the recent experimental possibility of Z=122 from 
natural 211,213 ' 217,218 Th-isotopes, associated with long lived 
superdiformed (SD) and/ or hyperdeformed (HD) isomeric 
states (UrOl], by 16 to 22 orders of magnitude longer than 
their corresponding ground-state, and more recently the syn- 
thesis of Z=117 at Flerov Laboratory 11611 from the reaction 
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20 Ca+ 2 l 4y Bk^ 29 '117 (Z=117, A=297), which decay simul- 
taneously via three and four neutrons into two differnet iso- 
topes 293 1 17 and 294 117, motivates us to focus on their prop- 
erties, using a microscopic theoritical model with better pre- 
dictive power. Such estimations of structure properties of nu- 
clei in the superheavy mass region is a challenging area in 
nuclear physics and a fruitful path towards the understanding 
of "Island of stability" beyond the spherically doubly-magic 
nucleus 208 Pb. 

The paper is organised as follows. Section II gives a brief 
description of the relativistic mean field formalism. The pair- 
ing effects for open shell nuclei, included in our calculations, 
are the same as discussed in |Q . The results of our calcula- 
tion are presented in Section III, and Section IV includes the 
a-decay modes of 293 1 17 and 294 1 17 isotopes. A summary of 
our results, together with the concluding remarks, are given in 
the last Section V. 



H. THE RELATIVISTIC MEAN-FIELD (RMF) METHOD 

The relativistic Lagrangian density for a nucleon-meson 
many -body system 117lll8ll . 
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All the quantities have their usual well known meanings. 
From the above Lagrangian we obtain the field equations for 
the nucleons and mesons. These equations are solved by ex- 
panding the upper and lower components of the Dirac spinors 
and the boson fields in an axially deformed harmonic oscil- 
lator basis, with an initial deformation (3q. The set of cou- 
pled equations is solved numerically by a self-consistent it- 
eration method. The centre-of-mass motion energy correc- 
tion is estimated by the usual harmonic oscillator formula 
E c .m. = |(41A -1 / 3 ). The quadrupole deformation param- 
eter /?2 is evaluated from the resulting proton and neutron 



quadrupole moments, as Q = Q n + Q p 



The root mean square (rms) matter radius is defined as (r^j = 
\ J p(r±, z)r 2 dr, where A is the mass number, and p(r±, z) 
is the deformed density. The total binding energy and other 
observables are also obtained by using the standard relations, 
given in lfl8ll . We use the well known NL3 parameter set lfl9tl . 
This set reproduces the properties of not only the stable nuclei 
but also well predicts for those far from the /^-stability valley. 
As outputs, we obtain different potentials, densities, single- 
particle energy levels, radii, deformations and the binding en- 
ergies. For a given nucleus, the maximum binding energy cor- 
responds to the ground state and other solutions are obtained 



as various excited intrinsic states. 

The constant gap, BCS-pairing approach is reasonably 
valid for nuclei in the valley of /3-stability line. However, this 
method breaks down when the coupling of the continuum be- 
comes important. In the present study, we deal with nuclei on 
or near the valley of stability line since the superheavy ele- 
ments, though very exotic in nature, lie on the /3-stability line. 
In order to take care of the pairing effects in the present study, 
we use the constant gap for proton and neutron, as given in 
ll20ll . which are valid for nuclei both on or away from the sta- 
bility line (for more details, see, e.g., Ref. lfl5ll . where E pa i r , 
the pairing energy, is also defined). 



III. RESULTS AND DISCUSSION 
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In many of our previous works and of other authors lUOl llE 
IT9I I2TI - I23I1 . the ground state properties, like the binding en- 
ergies (BE), pairing energies E pa i r , quadrupole deformation 
parameters f$i, charge radii (r c h), and other bulk properties, 
are evaluated by using the above stated relativistic Langragian 
for different forces. From these predictions, it is found that, 
generally speaking, most of the recent parameter sets repro- 
duce well the ground state properties, not only of stable nor- 
mal nuclei but also of exotic nuclei far away from the valley 
of /3-stability. This means to say that if one uses a reasonably 
well accepted parameter set, the predictions of the model will 
remain nearly force independent. In this paper we have used 
the successful NL3 parameter sets for our calculation. 
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FIG. 1: The binding energy per particle BE/A for the 284 -^ 10 H7 
isotopes, obtained in RMF(NL3) formalism and compared with the 
FRDM results 12411 . wherever available. 



A. Binding energy and two-neutron separation energy 

The calculated binding energy per nucleon BE/A and bind- 
ing energy BE, obtained from the RMF(NL3) formalism, are 



FIG. 2: The total binding energy BE for 284 - 310 H7 nuc i e i j n 
RMF(NL3) compared with the FRDM results |H. 
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FIG. 3: The two-neutron separation energy S2 n for 281 
clei, obtained from RMF(NL3) formalisms, and compared with the 
FRDM results 12411 . wherever available. 



comapred, respectively, in Figs. 1 and 2 and in Table I, with 
the Finite Range Droplet Model (FRDM) results J3. We 
notice that the BE/A obtained in the RMF(NL3) model over- 
estimate the FRDM result. In general, the BE/A value starts 
increasing with the increase of mass number A, reaching a 
peak value at A=288 for RMF(NL3) and at A=290 for the 
FRDM formalism. This means to say that 288 117 is the most 
stable isotope from the RMF(NL3) result and 290 117 from 
the FRDM predictions. Interestingly, 288 117 (with N=171) 
and 290 117 (with N=173) are both closer to the predicted 
closed shell at N=172 than at N=184. Note that the isotopes 
300 ' 302 117, next to the magic number N=184, are also in- 
cluded in this study. For the total binding energy BE of the 
isotopic chain in Table I and Fig. 2, we notice that the mi- 
croscopic RMF binding energies agree well with the macro- 
microscopic FRDM calculations, their differences decreasing 
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TABLE I: The RMF(NL3) results for binding energy BE, two-neutron separation energy S2n, pairing energy E pa ir, the binding energy differ- 
ence AE between the ground- and first-exited state solutions, and the quadrupole deformation parameter /?2, compared with the corresponding 
Finite Range Droplet Model (FRDM) results 12411 . The energy is in MeV. 





RMF(NL3) Result 


FRDM Result 


Nucleus 


BE 




Epair 


AE 




BE 


S2n 


0a 


288 


2052.586 


14.836 


14.698 


0.333 


0.018 


2047.09 


15.16 


0.080 


290 


2066.138 


13.552 


14.274 


0.360 


0.017 


2061.65 


14.56 


0.080 


292 


2079.802 


13.664 


14.109 


0.096 


-0.017 


2075.72 


14.07 


0.072 


294 


2092.468 


12.775 


13.653 


0.031 


0.041 


2089.22 


13.50 


-0.087 


296 


2104.803 


12.335 


13.583 


0.104 


0.028 


2102.66 


13.45 


-0.035 


298 


2116.598 


11.691 


13.274 


0.389 


0.015 


2114.79 


12.13 


-0.008 


300 


2128.174 


11.576 


12.841 


0.970 


0.005 


2126.14 


11.34 


0.000 


302 


2138.662 


10.488 


12.623 


0.596 


0.004 


2136.25 


10.11 


0.000 


304 


2148.296 


9.634 


12.695 


0.012 


0.002 


2145.71 


9.46 


0.000 


306 


2157.726 


9.430 


12.348 


0.004 


0.030 


2154.84 


9.13 


0.000 


308 


2167.327 


9.601 


11.912 


0.304 


0.047 


2163.93 


9.09 


0.001 


310 


2176.656 


9.329 


11.538 


0.512 


0.051 


2172.61 


8.68 


0.000 



gradually towards the higher mass region (around A=298), 
and then beyond this mass number the two curves again show- 
ing a similar behavior. Note that 298 117 (with N=181) is in 
this case closer to N=184. 

The two-neutron separation energy S , 2„(N,Z)= BE(N,Z)- 
BE(N-2,Z) is also listed in Table I. From the table, we find 
that the microscopic RMF Si n values also agree well with the 
macro-microscopic FRDM calculations. This comparison of 
Si n f° r RMF with the FRDM result are further shown in Fig. 
3, which clearly shows that the two S2n values coincide re- 
markably well, except at masses A=290 and 296. Apparently, 
the Sin decrease gradually with increase of neutron number, 
except for the noticeable kinks at A=290 (with N=173) and 
A=300 (with N=183) in RMF, and at A=296 (with N=179) in 
FRDM. Interestingly, these neutron numbers for RMF(NL3) 
are close to the earlier predicted 181 41011 N=172 or 184 magic 
numbers. 



B. Shape co-existence 

We have also calculated, for the whole Z=117 isotopic 
chain, the solutions in both prolate and oblate deformed con- 
figurations. In many cases, we find low lying excited states. 
As a measure of the energy difference between the ground 
band-head and the first excited state, we have plotted in Fig. 4 
(a) the binding energy difference AE between the two solu- 
tions, noting that the maximum binding energy solution refers 
to the ground state (g.s.) and all other solutions to the intrinsic 
excited state (e.s.). From Fig. 4 (a), we notice that in RMF cal- 
culations, the energy difference AE is small for the whole re- 
gion of the considered isotopic series. This small difference in 
the binding energy for neutron-deficient isotopes is an indica- 
tion of the presence of shape co-existence. In other words, the 
two solutions in these nuclei are almost degenerate for a small 
difference of output in energy. For example, in 290 117, the 
two solutions for /?2=0.017 and 0.123 are completely degen- 
erate with the binding energies 2066.138 and 2065.778 MeV. 
This later result means to suggest that the ground state can be 
changed to the excited state, and vice-versa, by a small change 



in the input data, like the pairing strength, etc., in the calcu- 
lations. In any case, such a phenomenon is known to exist in 
many other regions of the periodic table [25]. 
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FIG. 4: (a) The energy difference between the ground-state and the 
first-excited state AE, and (b) the pairing energy E pa ir, for the rel- 
ativistic RMF(NL3) calculation of Z=l 17 isotopic chain. 



Pairing is important for open shell nuclei whose value, for 
a given nucleus, depends only marginally on quadrupole de- 
formation Pi. This means that for differing /^-values in a nu- 
cleus, the pairing energy E pa i r changes only marginally (by 
~5-6%). On the other hand, even if the fa values for two 
nuclei are same, the E pa i r '$ could be different from one an- 
other, depending on the filling of the nucleons. This result is 
illustrated in Fig. 4 (b) for the RMF(NL3) calculation, where 
Epair for both the g.s. and first-excited state (e.s.), refering 
to different fa -values, are plotted for the full isotopic chain. 
It is clear from Fig. 4(b) that E pa i r decreases with increase 
in mass number A, i.e., even if the fa values for two nuclei 
are the same, the E pa i r 's are different from one another. This 
change of E pa i r is ~15% in going from, say, A=286 to 310. 
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FIG. 5: Quadrupole deformation parameter obtained from relativistic 
mean field formalism RMF(NL3), Compared with the FRDM results 
1 24], wherever available. 



FIG. 6: The rms radii r m of matter distribution and charge radii r c h 

for 286-310 U7 nudei; 

using the relativistic mean field formalism 

RMF(NL3). 



C. Quadrupole deformation parameter 

The quadrupole deformation parameter 02, for both the 
ground and first excited states, are also determined within 
the RMF formalism. In some of the earlier RMF and 
Skyrme Hartree-Fock (SHF) calculations, it was shown that 
the quadrupole moment obtained from these theories repro- 
duce the experimental data pretty well ITol [r7r - [l9l 

mm- 

l29ll . The g.s. quadrupole deformation parameter P2 is plotted 
in Fig. 5 for RMF, and compared with the FRDM results 12411 . 
It is clear from this figure that the FRDM results differ from 
the RMF(NL3) results for some mass regions. For example, 
the g.s. oblate solution appear for the nucleus 292 1 17 in RMF 
but is a prolate solution in FRDM. A more careful inespection 
shows that the solutions for the whole isotopic chain are pro- 
late, except at A=292 for RMF and at A=294-298 for FRDM 
model. In other word, there is a shape change from prolate 
to oblate at A=292 for RMF and at A=294 for FRDM. Inter- 
estingly, most of the isotopes are almost spherical in their g.s. 
configurations. 
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FIG. 7: (a) The Q a -energy and (b) the half-life time T a for a-decay 
series of 293 117 nucleus, using the relativistic mean field formalism 
RMF(NL3), compared with the FRDM data 13], the extrapolated 
experimental data I32U33I1 and the experimental data IT^l . whereever 
available. 



D. Nuclear radii 

The root-mean-square matter radius (r m ) and charge radius 
(T c h) for the RMF(NL3) formalism are shown in Fig. 6. As 
expected, the matter distribution radius r m increases with in- 
crease of the neutron number. However, though the proton 
number Z= 1 1 7 is constant for the isotopic series, the r c h value 
also increases with neutron number. A detailed inspection of 
Fig. 6 shows that, in the RMF calculations, both the radii 
show the monotonic increase of radii till A=310, with a jump 
to a lower value at A=292 (with N=175). There is no data or 
other calculation available for comparisons. 



IV. THE Q a ENERGY AND THE DECAY HALF-LIFE T a 

The Q a energy is obtained from the relation l30ll : 

Q a (N, Z) = BE(N, Z) - BE(N - 2, Z - 2) - BE(2, 2). 

Here, BE(N, Z) is the binding energy of the parent nucleus 
with neutron number N and proton number Z, BE(2, 2) is 
the binding energy of the a-particle ( 4 He), i.e., 28.296 MeV, 
and BE(N — 2, Z — 2) is the binding energy of the daughter 
nucleus after the emission of an a-particle. 
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TABLE II: The Q a energy and T a for Q-decay series of 117 nucleus, calculated on the RMF(NL3) model, and compared with the Finite 
Range Droplet Model (FRDM) results [24], the extrapolated data [32, 33] and the experimental data fl^l . wherever available. The energy is in 
MeV and the half-life time in second. 







RMF (NL3) Result 


FRDM Results 


Extrapolated result 


Experimental Results 


Nucleus 


Z 


BE 


Qa 


T 

± a 


BE 


Qa 


T 

± a 


BE 


Qa 


T a 


Qa 


T 

± a 


293 


117 


2086.602 


11.480 


-2.71 


2083.06 


11.68 


-2.40 








11.03 


-1.60 


289 


115 


2069.786 


10.552 


-0.96 


2066.45 


10.03 


1.26 


2063.17 


10.57 


-1.01 


10.31 


-0.31 


285 


113 


2052.042 


9.765 


0.60 


2048.18 


8.97 


3.86 


2045.45 


9.99 


-0.04 


9.74 


0.67 


281 


111 


2033.511 


9.231 


1.55 


2028.85 


9.37 


1.90 


2027.13 


9.69 


0.16 


9.48 


0.78 



TABLE III: Same as for Tabe II, but for 294 1 17 nucleus. 







RMF (NL3) Result 


FRDM Results 


Extrapolated result 


Experimental Results 


Nucleus 


Z 


BE 


Qa 


Ta 


BE 


Qa 




BE 


Qa 


Ta 


Qa 


Ta 


294 


117 


2092.578 


10.763 


-0.906 


2089.22 


11.66 


-3.13 








10.81 


-1.03 


290 


115 


2075.045 


11.083 


-2.325 


2072.59 


9.93 


0.77 


2069.730 


10.329 


-0.358 


9.95 


0.714 


286 


113 


2057.832 


9.612 


1.055 


2054.22 


8.90 


3.32 


2051.764 


9.752 


0.639 


9.63 


1.002 


282 


111 


2039.148 


9.080 


2.027 


2034.82 


8.80 


2.95 


2033.220 


9.464 


0.830 


9.00 


2.286 


278 


109 


2019.932 


8.775 


2.316 


2015.32 


9.38 


0.41 


2014.388 


9.176 


1.028 


9.55 


-0.099 


274 


107 


2000.411 


8.105 


3.913 


1996.40 


8.71 


1.81 


1995.268 


8.348 


3.042 


8.80 


1.517 
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FIG. 8: Same as for Fig. 7, but for 294 1 17 nucleus. 



of 117 are compared in Table II with the finite range droplet 
model (FRDM) calculation JH], the extrapolated lUllH] and 
the experimental data fl6Tl . wherever possible. The same com- 
parision is also carried out in Figs. 7 (a) and 7 (b), respec- 
tively, for Q a energy and the half-life time logioT a (s). 

From Figs. 7(a) and (b), and Table II, we notice that the 
calculated values for both Q a and T a (s) agree well with the 
known extrapolated as well as experimental data, but are over- 
estimated with-respect-to the FRDM predictions. For exam- 
ple, the values of T a from RMF coincide well for the whole 
mass region with the available experimental data, and with ex- 
trapolated values for 281 Rg, 285 113 and 289 115 nuclei. Sim- 
ilarly, for 2S1 Rg, the FRDM predictions match both the ex- 
trapolated and experimental results. Furthermore, the possible 
shell structure effects in Q a , as well as in T a (s), are notice- 
able for the daughter nucleus 285 113 (with N=172) for both 
the RMF predictions and experimental data. Note that N=172 
refers to the predicted magic number. 



With Q a energy at hand, we esimate the half-life 
logxoT a (s) by using the phenomenological formula of Viola 
and Seaborg Bltl: 



log w T a (s) = 



aZ — b 



(cZ + d). 



Here, Z is the atomic number of parent nucleus, and 
o=l. 66175, 6=8.5166, c=0.20228 and <i=33.9069. 



A. The a-decay series of 293 117 nucleus 

The binding energies of the parent and daughter nuclei are 
obtained by using the RMF formalism. From these BE, we 
evalute the Q a energy and the half-life time logioT a (s), using 
the above formulae. Our predicted results for the decay chain 



B. The a-decay series of 294 117 nucleus 

In this subsection, we present the Q a and the logioT a (s) 
results for decay series of 294 1 17 nucleus, using the same pro- 
cedure as in the previous subsection for 293 117. The results 
obtained are listed in Table III and plotted in Figs. 8(a) and 
8(b), compared with the FRDM predictions lf24ll . the extrapo- 
lated 1I32I, l33ll and experimental data |[l6ll . wherever possible. 

From Fig. 8(a) and (b), and Table III, we found almost 
similar results as are predicted in the previous subsection for 
293 117. Thus, the RMF(NL3) results for both Q a and T a (s) 
agree well with the known extrapolated and experimental data, 
but once again over-estimate the FRDM results. For example, 
the T a values for RMF coincide well with the experimental 
data for the whole isotopic chain and with the extrapolated 
data for 290 1 15, 286 1 13, 282 Rg, 278 Mt and 274 Bh nuclei. Sim- 
ilarly, FRDM predictions for 278 Mt and 290 1 15 match the ex- 
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trapolated and for 7 Bh and 78 Mt with experimental results. 
The possible shell structure effects in Q a , as well as in T a (s), 
are noticed for the daughter nucleus 286 113 (with N=173) for 
RMF and 278 Mt (with N=169) in experimental data, again co- 
inciding with earlier predicted N=172 magic number. 

V. SUMMARY 

Summarizing, we have calculated the binding energy, the 
rms charge and matter radii, and quadrupole deformation pa- 
rameter for the isotopic chain of recently synthesized Z=l 17 
superheavy element for both the ground- as well as intrinsic 
first-excited states, using the RMF formalism. From the calcu- 
lated binding energy, we have also estimated the two-neutron 
separation energy and the energy diffence between ground- 
and first-excited state for studying the shape co-existence, for 
the isotopic chain. Also, we have estimated the pairing en- 
ergy for the ground-state solution in the whole isotopic chain. 



We found a shape change from oblate to prolate deformation, 
with increase of isotopic mass number, at A=292. Most of the 
ground-state structures are with spherical solutions, in agree- 
ment with the FRDM calculations. From the binding energy 
analysis, we found that the most stable isotope in the series is 
288 117, which is close to predicted magic number at N=172. 
Our predicted a-decay energy Q a and half-life time T a match 
nicely with the available extrapolated and experimental data. 
Some shell structure is also observed in the calculated quanti- 
ties at N=172 and/ or 184 from RMF calculations of the vari- 
ous isotopes of Z=l 17 nucleus. 
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